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INTRODUCTION 

This study was undertaken to determine mission profile and delta 
velocity requirements to place a telescope at the Earth Su 
libration point L2 . The program "Integrated Mission Program 
(IMP), was selected to be used m the investigation. 

A description of IMP and its capabilities maybe found in 
Addenda ? The Addenda also contains the libration halo equ , 

constants and other parameters. Comments regarding the chaotic 
nature of numerical integration near the libration points are 
also attached in the Addenda. 

A basic two stage S/C with a simple mission profile was selected. 
This profile is shown immediately below. 

SPACE SURVEILANCE TELESCOPE MISSION PROFILE 
TBD DAY TRANSFER WITH COURSE CORRECTION AT 2-3 DAYS 

DEPART FROM NODE OF PARKING ORBIT (TLIB) 

160 NM : INCLINATION 23.45 DEG 

DROP 1ST STAGE 


MIDCOURSE 

HALO ORBIT INSERTION 
DROP 2ND STAGE 


The original study began with a s P» ce = r f £t ”® i 9 SL°two 

6200 lbs. This was raised scaling by 10 to 02000 Ips. The 
S/C are described in Table 3 and 6 respectively . y 

results will be shown for both vehicles. 

An important part of the study was to identify opportunites 
for departure and arrival. The initial study focused on a 
constant date of arrival. Trip times of 20 30 40, 60 80, 

and 100 days were investigated to arrive on 25 Sep 2011. 
fater a constant trip time of 100 days to arrive at selected 
dates in September, October, November, and December was studied. 

TRANSFER (TLIB) 

In a classical Newtonian two body system, the optimimum transfer 
from one circular orbit to another is a Hohmann transfer of 
degrees In an oblate system, an exact 180 degree transfer in 
inertial space cannot be accomplished by a single burn at the 
departure point, but must have a midcourse to account for the 
transfer plane's regression. Or a transfer can be made to the 
repressed apsides. Our system is a three body oblate system, 
wit (/perturbations that also pose restrictions .Two techniques 
for surmounting these restrictions are discussed on the next page 
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DEPARTURE -ARRIVAL POINTS 

Another requirement was an algorithm to ^ e S lfter in9 

departure date, time, orbit, position etc to arrive aft 
varying trip times. The departure and arrival points, Table , , 

were selected using the following as an algorithm. 

Select a transfer time and an arrival date. This choses a 
departure date and time. The departure 

program time zero, and sets the coordinate system that will 
be used for the calculations. 

With program zero time set, use the insertion event in 1MP 
to insert a target vehicle at the arrival time at the L2 point. 
From the output 9 note the arrival geocentric Latitude and the 
inertial Longitude. Note also whether the arrival azimuth is 
ascending or descending. The departure point is set to the 
opposite side of the Earth in plane with the Earth-Sun. 

Using IMP INSERTION, option 0, set the main S/C at the 
negative of the arrival Latitude, at the arrival inertial 
Longitude + 180 degrees, and at the departure altitude of 160 NM . 
Set 9 the azimuth flag to compute the proper direction of the 
orbit velocity. Input an inclination of 23.45 degrees. Then i 
velocity is set to 0 . , IMP will set circular velocity and azimuth 

This yields an almost inplane departure with close to 180 
decrees of transfer. A small coast can then set the transfer 
less than 179 degrees so that it can be solved. This coast 
also be used for optimization. A second technique to 
get ally ?rom the 180 Lgree point is to offset the longttude 
by a few degrees. This accounts for the transfer plane s 
regression rate varying as the S/C moves from the oblate 
gravity field to a spherical field. The second method was 
not used other than a few selected test cases. 

Table 1 assumes a constant arrival date 25 Sep 2011 is desired 
and shows departure Lat/Lon for several trip times. 

TABLE 1. CONSTANT ARRIVAL DATE 25 SEP 2011 E/S L2 


TRIP TIME 


DAY 

20 

30 

40 

60 

80 

100 


DATE 
SEP 05 
AUG 
AUG 
JUL 
JUL 
JUN 


26 

16 

27 

07 

17 


DEPARTURE 

LAT LON (G) 

- .731 107.68 DSC 
» 117.54 DSC 

» 127.40 DSC 

« 147.11 DSC 

" 166.82 DSC 

« 186.53 DSC 


ARRIVAL 


DATE 
SEP 25 


LAT 

.731 


LON (IN) 
-72.32 
-62.46 
-52.60 
-32 .89 


ASC 


-13 

6 


18 

53 
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Table 2 assumes a constant trip time of 100 days to arrive 
at the dates shown. 


TABLE 2 

DATE 
JUN 7 
JUN 17 
JUL 07 
AUG 07 
SEP 07 


100 DAY TRIP MIXED ARRIVAL DATES E/S L2 


DEPARTURE 
LAT 
3 . 146 
- .731 
-8.398 
-18 .399 
-23.298 


187 

186 

185 

184 

187 


LON (G) 
.418 DSC 
.530 DSC 
.043 DSC 
.691 DSC 
.400 DSC 


DATE 
SEP 15 
SEP 25 
OCT 15 
NOV 15 
DEC 16 


ARRIVAL 
LAT 
-3.146 
.731 
8.398 
18 .399 
23.298 


LON (IN) 
7.418 
6.532 
5 . 043 
4.691 
7 .400 


TLIB TRANSFERS TO L-2 

A simple method of breaking the initial orbit and moving 
out of the Earth's major influence was desired. Two 
stages are used. It is planned to use a fixed ■ 
stage for the TLIB burn, drop it and then after 2_3 days 
do a midcourse correction with the second stage. T e 
are several techniques for the TLIB burn. 

Departure point is perigee. The target apogee altitude 
is near the L2 radius from the Earth and. can be used 
for optimization. Velocity at departure is calculated 
from closed form two body equations then modif ^? T b 
BELL LAB'S correction for an oblate Earth. T 
burn is in the initial departure orbit plane. 

An inplane burn is made such that the S/C is in 
an orbit intersecting the departure orbit. At th 
intersection, the new orbit has a velocity and 
fliqht path angle as input. Normally 10908 M/S 
and 0 degrees. Again the begin event time is as 
given above and can be used for optimization. 


LAMBERT 

The halo inserti 
departure point 
Lambert problem 
initial point to 
desired. This so 
cannot be a 180 
may be modified 
algorithm "NXDN" 
This is referred 
the oblate Earth 
Since a Lambert 
the begin event . 
transfer. It sti 


on point is used as a target . The 
is the initial position. A spherical 
is solved for the velocity needed at the 
arrive at the target point at the time 
lution is basically for two bodies and 
degree inertial transfer. The transfer 
for the presence of the sun using IMP'S 
which numerically solves for a new XDN. 
to as an oblate Lambert . It includes 
gravity effects, and also the Sun's, 
solution for 180 degrees does not exist, 
state must insure less than 180 degrees 
11 may be used for optimization. 


OTHER THRUST EVENTS 
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MIDCOURSE . , . . 

An oblate Lambert is solved from the midcourse _ point _ 
to the halo insertion point to arrive at the time desire 
The midcourse transfer angle should be between 5 to 20 
degrees. And in a posigrade direction. Midcourse shou 
be accomplished while the Earth's gravity field still 
predominates. If delayed, things might get squirrely 
in the simulation. 

INJECTION , , , , 

Halo velocity calculated, and burn guided such that 
a minimum state error is obtained. On arrival this 
might require a big delta V. It is very affected 
by the length of time the sun's gravity has to 
change the approach trajectory. Injection is into 
a retrograde Halo at a point on the L0S . fr °™ ^2 
to the Earth. A current initial halo x is 2000 KM. 


NUMERICAL RESULTS 

The Integrated Mission program described on page 11 was used 
to integrate the equations of motion. IMP uses double precis 
FORTRAN and early in the study a requirement for at least 
Quad precision was suspected. Since this was not available, 
means of getting better results were investigated. 

See paqes 13 to 19 for a description of the basic equations 
used in IMP'S libration simulations. In particular note t e 
definition of a libration point. At such a point, the sum of 
an forces acting on a bod? is zero. Additionally attraction 

of gravity by either the sun or the earth is very small. Thi 
means that numerical roundoff and truncation can be a chao 
factor, and indeed, they were found so. 

One method for easing this problem is a modified patch conic 
procedure for TLIB and midcourse. The transfers are integrated 
with the Sun's qravity omitted from the equation of motion. 

Sn ?ea?Wng the 9 halo Insertion point, it is turned back on to 
calculate the Halo velocity needed. In the tables below, wit 
Sunref ers to results leaving the Sun in the equation of motion. 
A notation of OK means that the transfer did not exceed the 
capabilty of either stage. Most of the light vehicle transfe 
were in the not OK status. It was expected to be difficult to 
design a small system to put a telescope at L2 . 


TABLE 3 . GENERIC LIGHT VEHICLE DATA 

1ST STAGE: BREAK ORBIT AND TRANSFER 
THRUST 24000 LB 

ISP 300 SEC 

FUEL 4100 LB 

DRY 825 LB 

TOTAL 
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3185 M/S DV 
80 LB/S 


4925 LB 


2ND STAGE: MIDCOURSE AND FINAL INJECTION 528 M/S DV 

- ^ ~ ~ T D 


THRUST 

ISP 

FUEL 

DRY 

TOTAL 


1000 

285 


220 

325 


LB 

SEC 

LB 

LB 


3 .508 LB/S 


F45 LB 


PAYLOAD 

TELESCOPE 

CONTINGENCY 

TOTAL 

TOTAL INITIAL WEIGHT 


Basic HOHMANN transfers to radii near L2 were used to size 
the first stage of both vehicles. 

TABLE 4 . HOHMANN TRANSFERS LIGHT VEHICLE 

TRANSFER 
DAYS 
37 . 8 
41.7 

45.6 OK 

49.7 
53 . 8 



TABLE 5 . 

HOHMANN 

TRANSFERS 

LIGHT VEHICLE 


HOHMANN 

M/S 

WITH SUN G 

TDAYS 

TLIB 

MID 

INS 

SUM 

20 

3178 

98 

607 

3883 

30 

3178 

74 

378 

3629 OK 

40 

3178 

75 

359 

3612 OK 

60 

3176 

192 

384 

3752 

80 

3178 

329 

612 

4120 

100 

3178 

398 

715 

4292 


RADIUS 

PERIOD 

KM 

HRS 

1500000 . 

1818.4 

1600000 . 

2001.6 

1700000 . 

2190.5 

1800000 . 

2385 . 1 

1900000 . 

2585 . 0 


DV 

FUEL 

M/S 

LB 

3174.7 

4092.6 

3176.2 

4093.7 

3177 . 5 

4094.6 

3178 . 7 

4095.5 

3179 . 8 

4096.2 


700 

30 


LB 

LB 


73 0 LB 
6200 LB 
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TABLE 6. LIGHT VEHICLE OTHER TRANSFERS 


INPLANE 1ST FIX] 
TDAYS TLIB MID 

20 3175 104 

40 3175 106 

60 3175 220 

80 3175 366 

100 3175 399 

INPLANE 1ST FIX 
TDAYS TLIB MID 

20 3175 72 

40 3175 6 

60 3175 16 

80 3175 40 

100 3175 58 


M/S 

WITH SUN 

G 

INS 

SUM 


609 

3888 


329 

3610 

OK 

405 

3800 


574 

4115 


701 

4275 


M/S 

WITHOUT 

SUN G 

INS 

SUM 


548 

3795 


254 

3435 

OK 

330 

3522 

OK 

399 

3614 

OK 

451 

3685 

OK 


LAMBERT SPH M/S 


TDAYS . 

TLIB 

MID 

INS 

20 

3202 

20 

611 

30 

3237 

40 

369 

40 

3291 

81 

320 

60 

3345 

246 

426 

80 

3375 

427 

592 

100 

3402 

142 

493 


LAMBERT 

SPH 

M/S 

TDAYS 

TLIB 

MID 

INS 

20 

3202 

27 

552 

30 

3237 

28 

295 

40 

3290 

27 

254 

60 

3334 

27 

330 

80 

3378 

27 

402 

100 

3401 

26 

453 


WITH SUN G 
SUM 
3834 
3646 
3692 
4017 
4197 
4037 

WITHOUT SUN G 
SUM 
3781 
3559 
3571 
3702 
3807 
3880 
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TABLE 7. GENERIC HEAVY VEHICLE DATA 

1ST STAGE: BREAK ORBIT AND TRANSFER 
THRUST 60000 
ISP 300 

FUEL 41000 

DRY 3250 

TOTAL 

2ND STAGE : MIDCOURSE AND FINAL INJECTION 
THRUST 10000 
ISP 285 

FUEL 4250 

DRY 1250 

TOTAL 


'SFER 

3185 

M/S DV 

LB 

. 968 

T/W 

SEC 

200 

LB/S 

LB 



LB 




49250 LB 


[JECTION 

1130 

M/S DV 

LB 

.784 

T/W 

SEC 

35 

LB/S 

LB 



LB 




5500 LB 


PAYLOAD 


TELESCOPE 

CONTIGENCY 

TOTAL 


7000 

250 


LB 

LB 


PF .117 


TOTAL INITIAL WEIGHT 

See Tables 1 and 2 for the departure and arrival positions 
that were used for the following transfers. 

TABLE 8 . HEAVY VEHICLE TRANSFERS FIXED ARRIVAL DATE 


7250 LB 
62000 LB 


TDAYS 

20 

40 

60 

80 

100 


INPLANE TLIB 


TLIB 

MID 

3184 

45 

3184 

105 

3184 

231 

3184 

388 

3184 

401 


ARRIVAL 

SEP 25 

INS 

SUM 

614 

3843 

319 

3608 

433 

3848 

600 

4172 

715 

4301 


M/S 

WITH SUN G OK 
WITH SUN G OK 
WITH SUN G OK 
WITH SUN G OK 
WITH SUN G OK 


20 3184 

40 3184 

60 3184 

80 3184 

100 3184 


74 
6 

29 

55 

75 


552 

252 

336 

408 

459 


3810 

NO 

SUN 

G 

OK 

3492 

NO 

SUN 

G 

OK 

3549 

NO 

SUN 

G 

OK 

3648 

NO 

SUN 

G 

OK 

3719 

NO 

SUN 

G 

OK 


HEAVY INPLAE 



CD ZD IE 


s:\Ln x 'r-i o x m 


TRF DAY X 10**2 





TABLE 9. HEAVY VEHICLE 100 DAY TRANSFERS 


ARRIVE 

SEP 15 
SEP 25 
OCT 15 
NOV 15 
DEC 16 

SEP 15 
SEP 25 
OCT 15 
NOV 15 
DEC 1 6 
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INPLANE TL IB 


TLIB 

MID 

M/S 

3184 

273 

3184 

284 

3184 

268 

3185 

261 

3185 

257 

3184 

58 

3184 

58 

3184 

56 

3185 

50 

3185 

45 


INS 

SUM 

727 

4183 

727 

4195 

734 

4186 

742 

4188 

749 

4191 

449 

3691 

450 

3693 

463 

3703 

460 

3695 

465 

3695 


WITH SUN G 
OK 
OK 
OK 
OK 
OK 

NO SUN G 
OK 
OK 
OK 
OK 
OK 


HEAVY INPLANE 100 DAY 1=SUN 2=N0 



ARR MO X 10x*2 
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ADDENDA 

* * * INTEGRATED MISSION PROGRAM * * * 
* * * IMP * * * 

DEVELOPED AND PROGRAMMED 
BY 

V. A. DAURO, SR. 


* * * ABSTRACT * * * 

., IM p, IS A SIMULATION LANGUAGE THAT IS USED TO MODEL MOST PRESENT^ 

£ Si ss 

OR ?HR?E s?acSf? MAY BE USED ; A MAIN, A TARGET AND AN OBSERVER. 

* a FFHTRFRG 7/13 RUNGE-KUTTA INTEGRATOR WITH ERROR AND STEP SIZE 
CONTROL^ IS USED toUeRICALLY INTEGRATE THE EQUATIONS OF MOTION. 

* nnTATF DR FPHERICAL GRAVITY CAN BE USED FOR THE CENTRAL BODY. 
ADDITIONAL EFFECTS OF SUN GRAVITY, SOLAR PRESSURE, OR MOON GRAVIT1 
ARE IvAILABLE WHEN ADDED, THE SUN OR MOON GRAVITY IS SPHERICAL. 

* AERODYNAMIC LIFT AND DRAG WHILE IN THE EARTH OR MARTIAN ATMOSPHERE 
INCLUDED WHEN REQUESTED. 

- INPUT/OUTPUT HAS BEEN SIMPLIFIED AND IS IN MSTR^C ITOITS, WITH 
THE EXCEPTION OF THRUST AND WEIGHT WHICH ARE IN ENGLISH UNITS. 

™PU? 5s READ FROM THE VDT KEYBOARD AND THE USER'S INPUT FILE. 

REAL TIME KEYBOARD INPUT HAS BEEN MINIMIZED. 

* THE CODE EXECUTES EVENTS IN THE ORDER THEY ARE INPUT. AN EVEN* 

Tq FNTERED BY INVOKING ITS NAME, READING ITS OPTIONS, DATA, 
AND^THEN^EXECUTTNG . ON COMPLETION, THE NEXT EVENT IS STARTED. 

* EVENTS/MANUEVERS MAY INVOLVE NONE, ONE, °R ' 

tudt tqt t t AND PROPELLANT USAGE ARE AS THE USER SEL 

mraS'iTK IMPULSIVE, OR OF FINITE DURATION WITH 
™5 dm55 Slculated INTERNALLY OR PRESET BY THE USER. ^GORITHMS 
FOR TWO POINT BOUNDARY VALUES PROBLEMS INVOLVING VELOCITY CHANG ~ 

AND GUIDANCE ARE AUTOMATICALLY INVOKED AS NEEDED. 

„ MATM OUTPUT IS TO A USER NAMED PRINT FILE, TO A PLOT FILE, AND^TO 
A DEBUG FILE MAJOR MANUEVERS ALSO USE THEIR OWN PRENAMED FILEs 
?0 “^ ADDITIONAL DATA IN TABULAR AS WELL AS PLOT FORM. 

* THE CODE IS PROGRAMMED IN DOUBLE PRECISION FORTRAN AND WILL COMPI^a. 
WITHOUT^ ERRORS TO THE STANDARDS OF LAHEY (TM) FORTRAN 90. 

THE PROGRAM WAS INITIALLY CODED FOR MARSHALL SPACE FLIGHT CENTER , 
MEFU^ SE-AERO-G HUNTSVILLE, AL . THE AUTHOR WAS EMPLOYED BY NORTHRCX 

SERVICES rac HUNTSV I LLE , AL . LATER IT WAS EXTENSIVELY MODIFIED WHEN 

h“wM A MINION ANALYST AT PRELIMINARY DESIGN, PD33. MSFC. SINCE 

RETIREMENT HE HAS CONTINUED TO UPDATE AND IMPROVE IMP. 
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The following is data used and or computed in IMP to model 
the Lagrangian fibration point L2 . The spherical gravity. aode 
Of IMP and an ideal 3 body simulation ideal constants as below, 
were used to validate the equations of motion on page 15. 


* * * * EARTH / 

IDEAL LI BRAT I ON 
GM1 M**3/S**2 0 

GM2 M**3/S**2 0 

MEAN R M 
OMG SYS DEG/DAY 

DEG/S 0 

PERIOD DAYS 
OMG Ml DEG/SEC 0 


SUN * * * * 

YEAR 2011 DATA 

SYSTEM S/E 
398601200000E+15 
132718490000E+21 
149600000000 . 00 
0 . 98561057 
.114075297544E-04 
365.25582281 
, 417807419642E-02 


LI BRAT I ON 

RSYS KM 
OMG D/S 
RLIB KM 
VLIB M/S 
BC KM 


DATA SYS S/E 
YEAR MAX 
152098953 . 00 
0 . 00001170 
1526613 .13 
30338 . 524 
456.81 


POINT 2 

MIN 

147101055 . 15 
0 . 00001113 
1476449 .37 
29835 . 906 
441.80 


YEAR 2011 DATA 
DIF P/C 
4997897.85 3.29 

0.00000057 4.89 

50163.77 3.29 

502.618 1.66 
15.01 3.29 


PRESENT PARAMETERS AT LIB POINT 
GRAVITY PARTIAL WRT R 

EARTH -0.24767952E-12 

SUN - 0 . 8 0 924 3 8 9E - 13 


FORCE 
0 . 18284634E- 03 
0 . 60118772E-02 


(M/ S* * 2 ) 
(M/S * * 2 ) 


CENTRIFUGAL 


0 .61947235E-02 (M/S**2) 


DISTANCES 
SYSTEM 
LP TO EARTH 
BC TO SUN 


147103631.82 (KM) 
1476475.23 (KM) 
441.80 (KM) 
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EXTRACTED FROM IMP USER04 . MAN 
USER04 . MAN PAGES 14-22 

* * * LI BRAT I ON SIMULATIONS * * * 

LAGRANGIAN POINTS 

LIBRATION SIMULATIONS CAN BE CHAOTIC. CALCULATIONS, EVEN DOUBLE 
PRECISION, ARE NOT STABLE. ROUNDOFF AND TRUNCATION DEGRADE ACCURACY 
FOR S/E SYSTEM, OR E/M SYSTEM, IMP IS STABLE FOR A SIMULATION OF 
ABOUT ONE HALO ORBIT, S/E 180 DAYS, M/E 11 DAYS. USE LONGER TIMES 
WITH CAUTION. 

LIBRATION SIMULATIONS MAY BE 

1) AN IDEAL RESTRICTED THREE-BODY SOLUTION USING 
CIRCULAR ORBITS, SPHERICAL GRAVITY, IDEALIZED 
CONSTANTS, (MU, RADII, ETC) AND EITHER 

A) CLOSED FORM EQUATIONS OR, 

B) NUMERICAL INTEGRATION OF EQUATIONS OF MOTION 

2) A REAL WORLD SIMULATION USING 

ELLIPTICAL ORBITS, OBLATE GRAVITY, EPHEMERIDES , 

CURRENT CONSTANTS, (MU, RADII, ETC) AND EITHER 

A) APPROXIMATIONS USING PERTURBATION ANALYSIS OF 
THE IDEAL AND ASSUMED CLOSED FORM EQUATIONS 

B) NUMERICAL INTEGRATION OF EQUATIONS OF MOTION 

IT IS EVIDENT THAT THE ABOVE ARE ALL USEFUL TOOLS IN STUDYING 
LIBRATION. 1A IS USED FOR FAST PRELIMINARY RESULTS AND AT MOST 
IS ONLY ADVISORY. IB IS USED TO TEST THE IDEAL EQUATIONS AND THE 
RESULTS OF ROUND OFF AND TRUNCATION. 2A AGAIN GIVES FAST RESULTS 
AND IN GENERAL IS USED FOR PRELIMINARY PLANNING. 2B REQUIRES A 
LOT OF COMPUTER TIME AND IS NORMALLY USED LAST TO REFINE SOLUTIONS 
OBTAINED BY THE OTHER METHODS . 

IMP USES NUMERICAL INTEGRATION OF THE EQUATIONS OF MOTION, 

AND CAN BE USED IN AN IDEAL (IB) OR PERTURBED (2B) MODES. 
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KEATON'S EQUATIONS AS SHOWN IN REFERENCE 4, WERE ADAPTED AND 
USED IN IMP FOR LIBRATION STUDIES. REFERENCES 6, 7, AND 8 WERE 

ALSO USED IN PREPARING THE SIMULATION. THESE ARE ALL EXAMPLE^ 

OF THE CLOSED FORM IDEAL EQUATIONS. 

dv DEFINITION THE SUM OF FORCES AT A LIB POINT IS ZERO. EACH 
SYSTEM HAS 5 PLACES WERE THAT CAN OCCUR. THE FIRST THREE EARTH 
S UNBARE ON THE EARTH SUN LINE. LI ON THE SUNNY SIDE, U2 ON THE 
DARK SIDE, AND L3 ON THE SIDE OF THE SUN AWAY FROM THE EARTH 
ta tq ZXHFAD OF THE EARTH IN ITS ORBIT, AND L5 TRAILS THE EARTH 
IN 1?s“rb?t L™LS) SlTH THE EARTH AND SUN FORMS AN EQUILATERAL 
TRIANGLE. IMP HAS SIMULATIONS FOR LI, L2 , L4 AND L5 . L3 IS NOT 
NORMALLY STUDIED. 

IT IS KNOWN THAT THE LIB POINTS 1, AND 2 ARE UNSTABLE. HALO ORBITS 
aIoVT LI OR L2 CAN BE MADE STABLE IN THE RESTRICTED 3 BODY CASE. 

L4 AND L5 ARE STABLE AT WHAT ARE TERMED THE TROJAN POINTS. 

HALO ORBITS DO NOT EXIST ABOUT L4 AND L5 . 

HAT O ORBITS ABOUT LI AND L2 ARE USUALLY WHAT INVESTIGATORS STUD1 . 
WE INTEND STATIONING FOR A PERIOD OF TIME, IF PERTURBATIONS ARE 
PRESEN? INSTABILITY MAY RESULT. IN FACT, TRUNCATION OR ROUND OFF 
EVENTUALLY DISTURB THE SIMULATION OF THE SYSTEM. THEREFORE IN 
?ES?InG WE NEED TO EVALUATE THE IDEAL AS WELL AS THE REAL WORLD. 

TFSTS SHOW THAT IMP'S EQUATIONS, ARE ACCEPTABLE FOR AT LEAST 
ONE IDEAL EARTH/SUN, MOON/EARTH OR MARS/SUN HALO ORBIT. THAT IS, 

?N NON PERTTOBEi RESTRICTED 3 BODY MODE, THE PROGRAM WILL 

SIMULATE HALO ORBITS ABOUT E/S, E/M OR M/S LIBRATION POINTS LI 
AND L2 AS A MEASURE OF THE SIMULATIONS RELIABILITY, A HALO ORB I* 
ABOUT^ THE EARTH/ SUN LI POINT, RESTRICTED 3 BODY SOLUTION, NEEDS 
LESS THAN 1 M/S DELTA V TO ACCOUNT FOR ROUND OFF AND TRUNCATION. 

IDEAL STATIONING AT L4 AND L5 IS ALSO MODELED IN IMP, AND 90 
DAYS AT THE EARTH-SUN L5 POINT REQUIRES NO CORRECTIVE DELTA-V 

AS PART OF THE INVESTIGATION, PARAMETERS FOR THE IDEAL AND 
PERTURBED SYSTEMS WERE CALCULATED. FOR THE EARTH/SUN SYSTEM 
MODELED IN IMP , THE DATA OBTAINED WAS PREVIOUSLY SHOWN. 
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* * * PERTURBED results * * * 


STUDIES OF WHAT 
ECCENTRIC ORBITS 
SOME INTERESTING 


HAPPENS WHEN THE GRAVITY IS NON KEPLERIAN, OR 
ARE USED, OR A FOURTH BODY IS PRESENT YIELD 

RESULTS . 


w&T n ORRTTS ABOUT THE E/S LIBRATION POINTS LI AND L2 ARE FAIRLY 

Stable they are affected by the earth's orbit, and the presence of 

TOE MOON BUT THEY REQUIRE VERY LITTLE ENERGY TO REMAIN ON STATION. 

HALO ORBITS ABOUT THE M/E LIBRATION POINTS, ARE AFFECTED BY THE MOON ' S 
UNUSUAL ORBIT CONSIDER MOON EVECTION. WHENEVER THE MOON'S VELOTO. 
VfoTOEC?ED TOWARD THE SUN, IT SPEEDS UP. IT SLOWS DOWN WHEN MOVING 
> may "" FROM THE SUN . THIS EFFECT IS ALSO FELT BY THE 3RD BODY IN ITS 
HALO ORBIT THE SUN'S PRESENCE HAS AN EFFECT ON M/E HALO ORBITS ABOua 
ADUMBRATION POTO? STATIONING IS HOWEVER NOT AS DIFFICULT AS I ONCE 
THOUGHT. THE SUN'S EVECTION EFFECT IS EFFECTIVELY CANCELLED TO A 
GREAT DEGREE AND I NOW CONCLUDE THAT: 


* HAI O ORBITS ABOUT THE M/E LIBRATION POINTS (LI OR L2) CAN 
BE USED AS TRANSPORTATION NODES. DELTA-V FOR STATIONING 
IS APPROXIMATELY 7 TIMES THE EARTH SUN SYSTEM. 
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THE FOLLOWING EQUATIONS ARE MODELED IN IMP 

IDEAL LI BRAT I ON NODE/HALO SYSTEMS 
(SEE REFERENCES 4, 5, 6, 7, AND 8) 

B1 SMALLER BODY WITH GRAVITY GM1 

B2 LARGER " " GM2 

BC BARYCENTER OF SYSTEM 

LI LI BRAT I ON POINT B2 >>>> LI > B1 INSIDE 

L2 LI BRAT I ON POINT B2 >>>> B1 > L2 OUTSIDE 

DISTANCES 

RS B1 TO B2 SYSTEM 

R1 LI TO B1 

R2 LI TO B2 

R3 B2 TO BC 

R4 BC TO LI 

OMEGA SYSTEM ROTATION RATE 
OMEGA* * 2 = (GM1 + GM2 ) /RS* *3 


CALCULATIONS 
ALPHA = GM2/ (GM1+GM2 ) 

R3 = ( ONE- ALPHA) *RS 

ITERATE R1 UNTIL SUM OF F = VERY SMALL 


R2 = F (LIB POINT, Rl, RS) 
R4 = R2-RS 


Fl= GM 1 / R 1 * * 2 
F2= GM2 /R2 * * 2 
F3 = R4 *OMEGA* * 2 


GRAVITY BODY1 
GRAVITY BODY 2 

CENTRIFUGAL FORCE AROUND BC 


SIGN FI, F2 , AND F3 ACCORDING 
TO LIB POINT DESIRED AND SUM 


F1+F2+F3 = ZERO 


THEN SET KEATON'S EQUATIONS REFERENCE 4 
FSQ= (ONE-ALPHA) * (RS/R1) **3+ALPHA* (RS/R2 ) **3 
F=SQRT (FSQ) 

BSQ=ONE-FSQ/TWO+F*SQRT (2 . 25DO*FSQ-TWO) 

B=SQRT (BSQ) 

GAMH= (ONE + BSQ + TWO* FSQ ) / (TWO*B) 

BOMGS=B* OMEGA 

FOMGS=F* OMEGA _ 

SPER= 2*PI/0MEGA SYSTEM PERIOD 

HPER=SPER/B HALO PERIOD 
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COORDINATE SYSTEMS 

THE INTEGRATION OF THE DEQS OF MOTION IS DONE WITH THE 
BODY B1 AS THE CENTER OF THE COORDINATE SYSTEM. LI BRAT I ON 
NODES ARE SET IN THIS SYSTEM AND USED AS THE CENTER OF THE 
HALO COORDINATE SYSTEM (HX, HY , HZ) . 


HX POINTS FROM B1 TO B2 

HY IS TO LEFT VIEWED FROM TOP IN ROTATION PLANE 
HZ COMPLETES A RIGHT HAND SYSTEM 

TO INSERTION TIME (ONLY AT BANG = ZERO, OR PI) 

XO INSERTION POSITION ALWAYS POSITIVE 


YO ALWAYS AT ZERO 

ZO OUT OF PLANE COMPONENT AT 

PHASE 0 INSERTION PHASE ZERO 

PI 

T = POSITION TIME 

THALO = T - TO 
BANG=BOMGS*THALO+ PHASE 0 
FANG= FOMGS *THALO+ PHASEO 

HX= XO*COS (BANG) 

HY=-X0*GAMH* SIN (BANG) 

HZ= ZO*COS ( FANG) 

HXD=-XO*BOMGS*SIN (BANG) 

HYD= -XO*BOMGS*GAMH*COS (BANG) 
HZD=-Z0* FOMGS* SIN (FANG) 


INSERTION ALWAYS POSITIVE 
IF INSERT AT BANG = ZERO 
IF INSERT AT BANG = PI 


ANOMALY ANGLE IN HX, HY PLANE 
11 Z" ANGLE WRT TO HX, HY PLANE 
HALO POSITION 


HALO VELOCITY 
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* * * LI BRAT I ON NOTES * * * 


TT IS IMPORTANT THAT THE USER KNOW THE DIFFERENCE BETWEEN THE 
^qrllCTED THREE BODY SYSTEM (IDEAL) AND THE REAL WORLD IN 

THE IDEAL SYSTEM, THERE ARE ONLY THREE BODIES n^f^THE MASSES 
ARE IN IMP'S SIMULATION, DESIGNATED Bl, B2 , AND B3 . THE MASS w 

ARE SUCH THAT 


B3 < <B1<B2 . 

pi AND B2 ROTATE ABOUT A COMMON BARYCENTER WITH A CONSTANT RATE 

MINIMAL INFLUENCE ON THE SYSTEM. 

IN THE REAL WORLD, THERE ARE PERTURBATIONS. 

IN THE SUN-EARTH SYSTEM, THE EARTH IS IN ^_ E ^ P ^ C su C h R ?^T A THE T 
rap VPFMTER ALTHOUGH THE DISTANCE TO THE SUN IS SU 
If™ APPEARS^TO HAVE A SPHERICAL GRAVITY FIELD, THE THIRD BODY IS 
fff.SE ENoSgH TO THE EARTH SO THAT IT IS SEES AN OBLATE EARTH GRAVI 
FTFT.D FURTHER THERE IS A FOURTH BODY, THE MOON, CLOSE ENOUGH Tu 
flreCT TOE SYSTEM HOWEVER, FOR B3 AT THE EARTH-SUN LI BRAT I ON 
POINTS , ^THE MOON IS FAR ENOUGH AWAY SO THAT ITS EFFECTS ARE MINIMA 

FOR THE EARTH-MOON SYSTEM, THE SUN IS A FOURTH BODY THAT DOES 
AFFECT THE SYSTEM AND MUST BE ACCOUNTED FOR IN A SIMULATION. 
ADDITIONALLY THE MOONS ORBIT IS 3 TIME MORE ECCENTRIC THAN 
TOE EARTH, AND THE MOON ORBITS THE EARTH ABOUT 13 TIMES A YEAR. 
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